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SHEA^  STRESS-STRAIN  RELATION  OBTAINED  FROM  TORQUE-TWIST  DATA 
INTiiODUCTION 

One  method  of  establishing  the  shear  stress-strain 
relation  of  a  material  involves  applying  torsion  to  a 
circular  cylinder.  This  approach  has  the  advantage  of 
experimental  simplicity,  but  the  disadvantage  that  the 
specimen  is  nonuniformly  loaded,  thus  requiring  a  somewhat 
sophisticated  data  analysis.  Fortunately  the  required 
analysis  has  already  been  worked  out  for  the  case  of  solid 
cylinders  made  of  an  isotropic  ejAStic  material.  Such  an 
analysis  was  furnished  by  Nadai  (!].  Unfortunately,  some¬ 
times  only  an  approximate  analysis  is  performed,  and  results 
are  only  reported  as  the  resulting  approximate  shear  stress- 
strain  relation.  Fortunately  however,  as  we  shall  show,  this 
type  of  approximate  result  can  be  easily  corrected  to  get 
the  proper  result  without  having  to  go  back  to  the  original 
torque-twist  data,  which  may  no  longer  be  available. 

A  second  contribution  of  the  paper  is  to  extend  the 
analysis  to  the  case  of  hollow  cylinders.  A  third  is  to  call 
attention  to  the  fact  that  some  cases  of  combined  loads  can 
also  be  analyzed  in  essentially  the  same  manner.  This  is  of 
particular  interest  in  the  case  of  highly  anisotropic 
materials,  in  which  compression-shear  interaction  may  be 
large . 

We  first  derive  the  standard  result,  using  an  approach 
that  emphasizes  the  physical  aspects  of  the  problem,  and 
serves  as  the  basis  of  the  extension  to  hollow  cylinders. 

SOLID  ISOTROPIC  CYLINDER 

Consider  first  a  solid  circular  cylinder  in  uniform 
torsion,  in  which  the  twist  is  0  radians  per  unit  length. 
The  shear  strain  is  then  given  by 

y=r9  (1) 

where  r  is  the  radius  of  the  field  point  under  consi¬ 
deration.  If  the  cylinder  is  linearly  elastic,  the  shear 
stress  T  is  proportional  to  the  radius  as  shown  in  Fig.  1. 
If  the  material  yields  at  a  stress  below  the  maximum  shown 
in  the  figure,  a  plot  of  t  vs.  r  will  be  nonlinear.  Also 
included  in  the  figure  is  a  nonlinear  plot  for  the  same 
applied  torque. 

In  the  elastic  case,  the  analysis  is  quite  straight¬ 
forward.  The  torque  is  given  by 


M 


=  f  2nr  r  GrQJr  = 


na 


GQ 


(2) 


where  G  is  the  shear  modulus  of  the  material/  and  a  is 
the  radius  of  the  cylinder.  The  shear  stress  at  the  cylinder 
surface  is  given  by 

T(a  )  =  GQa  =  2M  Iva}  (3) 


An  approximate  shear  stress-strain  curve  for  the 
nonlinear  case  can  be  obtained  by  assuming  that  for  any 
value  of  twist  (3)  applies.  Since  (3)  is  only  valid  in  the 
elastic  case,  we  shall  refer  to  this  as  the  quasi-elastic 
(QE)  solution.  It  is  equivalent  to  replacing  the  curve  in 
Fig.  1  by  the  straight  line  and  finding  the  shear  stress- 
strain  relation  at  the  surface.  It  is  clear  that  this 
overestimates  the  shear  stress  for  a  given  strain,  but  it  is 
not  immediately  obvious  how  big  the  error  might  be. 

The  standard  result  [1-3]  is 


One  way  to  derive  this  result  is  as  follows.  In  the 
nonlinear  case,  for  twist  0  the  torque  is 


A/(9)  = 


=  |2jtr^(r 


where,  for  simplicity,  H  denotes  the  in+^.egrand.  When  the 
torque  is  increased  to  M  +  dM  the  twist  becomes  0  +  de  .  By 
virtue  of  (1),  every  stress  contour  moves  toward  the  center 
a  distance  equal  to  r  d0  /  9  .  This  is  depicted  graphically 
in  Fig. 2. 

It  is  convenient  to  consider  the  torque  for  M  ♦  dM  to 
be  composed  of  two  parts,  the  part  contributed  by  stresses 
below  T  (a,  9  )  and  the  part  contributed  by  stresses  larger 

than  this  value.  Thus 

a  ( 1  -  (^6(9)  a 

M(8-t-de)=  J  ndr  +  J  ndr  (6) 

1)  a(l  - 

Since  each  r  and  the  dr  in  the  first  integral  are  decreased 
by  the  factor  (1  -  de  /  e  )  while  the  product  r9  remains 

unaffected  by  the  increase  in  9  ,  we  can  write 

<1(1  -  <<e/8) 

/  nJr  =  (I  -  d9/9)^,W(e) 

0 

=  (7) 


Note  that  this  result  holds  regardless  of  the  form  of  t  ( y  ) • 
The  second  term  is  simply 


(8) 


a 

J  ndr  =  T(a  )2iw*d9/9 

a  (I  -  da/e) 

Combining  (5) -(8)  we  obtain 


(9) 

(4) 


From  (4)  it  is  evident  that  the  proper  reduction  of  the 
torque-twist  data  requires  knowing  both  the  ordinate  and  the 
slope  of  the  torque-twist  curve.  What  can  we  do  about  it  if 
we  only  have  the  QE  approximation  referred  to  earlier? 
Fortunately  the  correction  is  very  easily  made. 

The  approximate  result  was  obtained  by  assuming  that 
(3)  applies.  Let  us  designate  the  resulting  shear  stress  by 
the  symbol  .  Thus 

?  =  2W /jw^  (10) 


Then 


•c(a ) 


1 

2na^ 


Since  y 


=  t(a  )  -  -7  T(<j  ) 

4 

a9  we  finally  obtain 


(11) 


(12) 


It  is  readily  seen  from  Fig.  3  that  the  term  in  the  square 
bracket  in  (12)  is  simply  the  ^  -intercept  of  the  tangent  to 
the  approximate  (QE)  shear  stress-strain  curve  taken  at  any 
chosen  value  of  y  .  Below  the  elastic  limit  the  term  is 
zero.  The  term  increases  as  the  slope  decreases r  becoming 
equal  to  the  ordinate  when  the  slope  is  zero.  Assuming  the 
slope  never  becomes  negative  (at  least  for  values  of  y  that 
are  of  practical  interest)  the  correction  never  exceeds  25%. 
The  correct  curve  is  thus  bounded  by  the  two  curves  t  (a)  and 
0.75  ^ (a)  as  shown  schematically  in  Fig.  4. 
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HOLLOW  CYLINDER 


From  Eq.  (9)  we  see  that  in  the  case  of  a  solid 
cylinder  the  increment  in  torque  resulting  from  an  increment 
in  twist  is  composed  of  two  parts,  one  negative  and  one 
positive.  The  negative  one  is  associated  with  the  shrinking 
radius  of  each  of  the  shear  stress  contours,  while  the 
positive  one  is  the  consequence  of  a  higher  stress  being 
generated  at  the  outer  surface.  The  latter  contains  the 
surface  stress  t  (a, 9  )  as  a  factor.  On  physical  grounds  one 
can  anticipate  the  same  or  similar  terms  will  be  present  in 
the  case  of  a  hollow  cylinder;  additionally  there  will  be  a 
negative  term  containing  the  factor  t  (b,e )  because  when  the 
twist  increases  the  lowest  stress  zone  moves  out  of  the 
cylinder  from  the  inner  surface,  r  =  b. 

To  make  these  qualitative  considerations  quantitative, 
consider  a  hollow  cylinder  with  outer  and  inner  radii  equal 
to  a  and  b  respectively.  The  torque  M  required  to  twist  the 
cylinder  through  an  angle  9  per  unit  length  can  be  regarded 
as  the  difference  in  the  torques  (  9  )  and  M|j(9  )  required 
to  twist  solid  cylinders  of  radii^a  and  b  through  the  same 
angle.  That  is. 


M(Q)  =  A/a(0)  - 


(13) 


Now  from  (4)  the  exterior  surface  stresses  on  the  two 
solid  cylinders  are 
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(15) 

(16) 

(17) 


Comparing  (17)  with  (4)  it  is  evident  that,  as 
expected,  the  surface  stress  t  (a)  is  a  little  larger  for  a 
hollow  cylinder  than  for  a  solid  one  with  the  same  torque 
applied.  Also  as  b-^a,  t  (b)  -*■  t  (a)  ,  as>  expected. 

The  torque-twist  curve  of  the  holiuw  cylinder  enables 
us  to  evaluate  the  term  in  square  brackets.  But  this  yields 
only  a  single  equation  to  evaluate  the  two  unknown  stress¬ 
es  T  (a,  9  )  and  t  (b,  9  ) .  We  can,  nevertheless,  solve  for 
shear  stress  as  a  function  of  strain;  but  in  contrast  to  the 
solid  cylinder  case,  it  must  be  done  iteratively.  We  start 
by  choosing  a  value  of  9  such  that  T(b,0  )  is  known,  namely 
one  at  which  it  is  at  or  below  the  elastic  limit  x  .  In 
this  case,  ® 


I «. 


T(b,6)  =  Gbe  (18) 

If  G  is  not  already  known,  it  can  be  obtained  by  employing 
the  equation  found  by  changing  the  lower  limit  in  (2)  from  0 
to  b: 

M  =  j(a*  -  b*)G6  (19) 

If  the  next  angle  of  twist  is  chosen  to  be  larger  by  the 

factor  a/b,  the  corresponding  value  of  t  (b)  is  equal  to  the 
initial  value  of  t  (a) ,  i.e., 

1(6,02)  =  t(a, 01 )  (20) 

Repeating  this  process  one  will  always  know  the  value 
of  T  (b)  needed  to  evaluate  t  (a) ,  but  only  for  the  sequence 
of  twist  angles  6^^  =(a/b)  ^“^e^^.If  this  procedure  results  in 
the  calculated  points  on  the  shear  stress-strain  curve  being 
too  far  apart  and  more  points  are  desired,  we  can  proceed  as 
follows : 

Choose  0j^  to  be  equal  to  or  slightly  below  the  twist  at 
which  the  torque-twist  curve  becomes  nonlinear.  Then  choose 
successive  values  of  twist  to  satisfy  the  relation 

♦  1  =  0.  (a  lb  f'l  (21) 

where  j-1  is  the  number  of  points  desired  between  successive 
pairs  of  points  obtained  by  the  previous  procedure.  Note 
that  for  i  ^  j,  T(b,9j^)  is  at  or  below  the  elastic  limit 
and  thus  is  known  from  (18).  The  corresponding  values  of 
outer  surface  shear  stress  t  (a, 9  )  can  thus  be  found  from 

(17) .  For  i> j  , 

i(  6 ,0, )  =  T(  a  ,0i  )  (22) 

so  that  each  t  (b)  is  equal  to  a  previously  calcula¬ 
ted  T  (a)  . 

By  following  the  above  procedure,  t  (a  ,  e  )  can  be 
evaluated  over  the  entire  range  of  twist  9  .  The 
corresponding  shear  strain,  of  course,  is  equal  to  a9 

DISCUSSION 

The  discussion  up  to  this  point  has  assumed  that  the 
material  under  consideration  is  isotropic.  Actually,  of 
course,  the  procedure  is  equally  applicable  to  planar 
isotropic  materials,  provided  that  the  cylinder  axis  is 
normal  to  the  plane  of  isotropy.  The  shear  stress-strain 
curve  obtained  in  this  case  is  alternative¬ 
ly  ijy  7zy  where  the  x  and  y  axes  lie  in  the  plane  of 

isotropy. 

We  note  in  closing  that  although  the  preceding  analysis 
assumed  that  the  applied  axial  stress  is  zero,  it  continues 


to  hold  for  any  other  constant  value  of  axial  stress.  That 
is  because  in  this  case  the  combined  stress  contours  move  in 
toward  the  center  of  the  cylinder  in  just  the  same  way  as 
they  do  for  pure  torsion.  A  cylinder  torsion  test  under 
constant  applied  axial  load  thus  offers  a  good  experimental 
means  for  testing  shear-compression  coupling  in  planar- 
isotropic  materials.  This  coupling  is  currently  quite 
inadequately  understood;  but  it  can  be  very  important  when 
the  compressive  modulus  is  orders  of  magnitude  greater  than 
the  shear  modulus. 

CONCLUSIONS 

1.  A  new  derivation  has  been  given  for  the  correct  procedure 
for  converting  the  torque-twist  relation  into  a  shear 
stress-strain  relation  in  the  case  of  a  solid  isotropic 
cylinder.  While  the  result  is  not  new,  the  derivation  is 
based  more  on  physical  considerations  and  less  on 
mathematical  ones  than  most  previous  derivations,  and 
therefore  may  be  preferred  by  some  readers. 

2.  It  is  shown  that  results  sometimes  reported  using  the 
quasi-elastic  (QE)  approximation  for  this  purpose  can  be 
easily  processed  to  give  the  correct  result. 

3.  A  f)rocedure  is  described  for  similarly  finding  the  shear 
stress-strain  curve  from  the  the  torque-twist  curve  obtained 
usin']  a  hollow  cylinder. 

4.  The  procedure  described  is  also  applicable  to  the 
analysis  of  planar-isotropic  materials  provided  the  cylinder 
axis  is  normal  to  the  plane  of  isotropy. 

5.  Cylinder  torsion  tests  provide  a  potentially  valuable 
means  for  studying  shear-compression  interaction  in  highly 
anisotropic  planar-isotropic  materials. 
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